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INTEGRATING-MATRIX METHOD FOR DETERMINING 
THE NATURAL VIBRATION CHARACTERISTICS 
OF PROPELLER BLADES* 

By William F. Hunter 
Langley Research Center 

SUMMARY 

A numerical method is presented for determining the natural lateral vibration char- 
acteristics of a rotating twisted propeller blade having a nonuniform, unsymmetrical cross 
section and cantilever boundary conditions. Two coupled fourth-order differential equa- 
tions of motion are derived which govern the motion of such a beam having displacements 
in two directions. A development of the integrating matrix, which is the basis of the 
method of solution, is given. By expressing the equations of motion in matrix notation, 
utilizing the integrating matrix as an operator, and applying the boundary conditions, the 
differential equations are integrated and formulated into an eigenvalue problem whose 
solutions may be determined by various methods. Numerical examples are presented, 
and the computed results are compared with experimental data and exact solutions. 

INTRODUCTION 

The determination of the natural vibration characteristics is of fundamental impor- 
tance in the design of propeller blades. Propellers often have serious resonant vibration 
problems with the excitation frequencies usually being equal to the rotational speed or 
some multiple of it. To insure that conditions susceptible to resonance do not exist within 
the range of operating speeds, it is necessary that the natural frequencies be determined 
accurately as a function of the rotational speed. Also, the natural modes, because of 
their orthogonality relationships, are very suitable for use in "series” solutions of 
response problems. 

This paper formulates a numerical solution for the natural vibration frequencies 
and the modal functions of twisted propeller blades. In addition to presenting a means 

*Most of the information contained herein was included in a thesis entitled "Inte- 
grating Matrix Method for Determining the Natural Vibrations of a Rotating, Unsymmetri- 
cal Beam With Application to Twisted Propeller Blades" submitted in partial fulfillment of 
the requirements for the degree of Master of Science in Engineering Mechanics, Virginia 
Polytechnic Institute, Blacksburg, Virginia, June 1967. 
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for determining accurately the vibration characteristics of propeller blades, this paper 
presents a numerical method for solving linear ordinary differential equations. 

The twisted propeller is idealized in this analysis as a rotating cantilevered beam 
which has a nonuniform and unsymmetrical cross section. The term "twist" is used to 
define a variable orientation along the length of the beam of the principal axes relative to 
the plane of rotation and is not to be confused with torsion. Either twist or an unsym- 
metrical cross section causes the beam to have transverse displacements in two direc- 
tions. Two coupled differential equations, which describe the behavior of such a beam 
having motions in two directions, and the associated orthogonality relationships are 
derived in the appendixes. 

A number of studies have been published concerning the determination of the natural 
vibration characteristics of propeller blades by energy methods, by adaptations of the 
Holzer method, and by the Stodola method. Earlier analyses neglected the twist of the 
blades and also assumed the cross sections to be symmetrical. More recent studies 
(refs. 1 to 4) have treated the twisted propeller which has lateral deflections in two direc- 
tions. In reference 1 this problem was analyzed through the use of energy principles. 
Extensions of the Holzer method were used in references 2 and 3. Reference 4 gives a 
solution based upon the Stodola method. 

In this paper the integrating matrix is developed as the basis for the method of 
solution. The integrating matrix is a means of numerically integrating a function that is 
expressed in terms of the values of the function at equal increments of the independent 
variable. It is derived by expressing the integrand as a polynomial in the form of 
Newton's forward-difference interpolation formula. Integrating matrices based upon 
polynomials of degrees one to seven are given. 

The solution to the governing differential equations of motion is developed entirely 
in matrix notation. First, the differential equations, which are fourth-order linear homo- 
geneous equations having variable coefficients, are expressed in a matrix equation. This 
matrix differential equation is then integrated repeatedly by using the integrating matrix 
as an operator. Next, the constants of integration are evaluated by applying the boundary 
conditions. Finally, the resulting matrix equation is expressed in the familiar concise 
form of the eigenvalue problem. The solutions to this eigenvalue problem may be 
obtained by any one of several methods. Two methods which have been utilized for this 
problem are described. The solutions define the natural frequencies and the modal dis- 
placements. The dominant eigenvalues correspond to the lower frequencies. 

The continuous problem is treated instead of a discrete or lumped mass system. 
Integral equations as such and influence coefficients, which are typical of other integra- 
tion methods for boundary -value problems, are not used. In developing the solution, all 
functions are in effect represented by high-degree polynomials at the boundaries as well 
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as elsewhere within the beam. Since the polynomials approximate the continuous func- 
tions very accurately, the integration of these polynomial representations yield extremely 
small errors. In addition to treating the continuous problem and yielding accurate 
results, the method is appealing because the numerical solution may be formulated quickly 
from the differential equations and may be easily programed for computations by a digital 
computer. Also, the inputs are generated very simply since they are merely the coeffi- 
cients appearing in the differential equations. 

Numerical examples are presented to show that the solution is applicable to propel- 
ler blades and to give an indication of the accuracy of the method. The computed natural 
vibration frequencies of a rotating propeller blade are compared with experimental data. 
Also, the numerical solutions for the natural vibration characteristics of a nonrotating 
cantilevered beam, which has a uniform and symmetric cross section, are compared with 
the exact solutions. 


SYMBOLS 

A cross-sectional area, in^ (m2) 

a acceleration vector 

aj polynomial coefficient where i = 0, 1, 2, . . .,r 

a^ k element in the jth row and kth column of the matrix £aJ where 

j,k = 1, 2, . . ., n 

Cj constant of integration 

Cj scaling constant 

Di determinant of the matrix M with the ith and (n+i)th rows and the ith and 

(n+i)th columns removed 

E shifting operator; modulus of elasticity, lb/in2 (N/m2) 

f arbitrary function of x 

g function of p 

h increment of the independent variable x; length of beam interval, in. (m) 
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moments and product of inertia, in^ (m^) 


i station number 

i, J, k unit vectors along the X-, Y-, and Z-axes, respectively 

j index denoting jth natural mode of vibration 

k^,kjj,kj 2 constants of integration 
l length of beam, in. (m) 

M bending moment, in-lb (m-N) 

^ minor of the element aj ^ 

m mass per unit length, Ib-sec2/in2 (N-sec2/m2) 

n number of intervals 

p nondimensionalized independent variable 

R radius of propeller blade, in. (m) 

r degree of polynomial approximation; index denoting the rth natural mode 

of vibration 

T position vector 

s index denoting the sth natural mode of vibration 

T axial tensile force, lb (N) 

t time, sec 

u elastic displacement of the centroidal axis in the X-direction, in. (m) 

u total elastic displacement in the X-direction of a point not on the centroidal 

axis, in. (m) 
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V shear force, lb (N) 

v,w components of the elastic displacement of the centroidal axis in the direc- 

tions parallel and normal, respectively, to the plane of rotation, in. (m) 

v,w vibration amplitude for v and w, respectively, in. (m) 

w variable in numerical example 

X,Y,Z axes of rotating orthogonal reference frame 

x independent variable; longitudinal body-axis coordinate, in. (m) 

y,z cross-sectional body-axis coordinates which are parallel and normal, 

respectively, to the plane of rotation when the propeller is in the unstrained 
condition, in. (m) 

y’ ,z' inclined cross-sectional body-axis coordinates, in. (m) 

Q! j ^ cofactor of the element ^ 

/3 orientation of the y' ,z* reference with respect to the y,z reference, 

deg or rad 

A^Iq kth forward difference 

e longitudinal strain, in./in. (m/m) 

? arbitrary function 

A. eigenvalue, sec^ 

S2 angular velocity of the propeller blade, rad/sec or rpm 

w natural vibration frequency, rad/sec or Hz 

w angular velocity vector 
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Subscripts: 


i denotes function evaluated at x = Xj where i = 0, 1,2, . . .,n 

j denotes jth natural vibration mode 

r denotes degree of polynomial approximation 

Y,Z denote directions of moments and shear forces 

y,z,y* ,z' denote axes about which moments and products of inertia are taken 

Matrix notation: 


[] 

square matrix 

n 

diagonal matrix 

o 

column matrix 

LJ 

row matrix 

[J 1 

inverted matrix 

1 1 

' H 

transposed matrix 

HIT] 

integrating matrices 

M 

unit or identity matrix 

© 

column matrix with all elements unity 

(o) 

column matrix with all elements zero 


Matrices and the equations by which they are defined are as follows: 
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Matrix 

Equation number 

M 

(C7) 

M 

(Al) to (A7) 

[b] 

(10c) 

L B °J 

(F9a) 

H 

(F9b) 

M 

(22e) 

N 

(34a) 

W 

(34b) 

H 

(41) 

[f] 

(FI 3b) 

CO 

(8c) 

[°] 

(45c) 

[H] 

(45d) 

CO 

(12) 

[ 1 

(24a) 

( K )> 

(27), (35), (36), (37), (38) 

M 

(C4c) 

H 

(22b) 


(F3d) 
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Matrix 


[n] 

[n(s)) 

Kl 

L Q J 

DO 

DO 

W 

M 

{t} 

w 

H 

CO 

CO 


OD 




Equation number 
(46b) 

(E8) 

(43b) 

(Elb) 

(2 2d) 

(F3c) 

(C4b) 

(22c) 

(39b) 

(F3b) 

(28b) 

(E2) 

(22a) 

(E6a) 

(8b) 

(10b) 
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ANALYSIS 


In the analysis, a development of the integrating matrix is presented. A statement 
of the governing differential equations and boundary conditions is given. The differential 
equations are expressed in matrix notation and are then integrated numerically by using 
the integrating matrix as an operator. After the evaluation of the constants of integra- 
tion, an eigenvalue problem results. Also, procedures for solving the eigenvalue problem 
are described. 


Development of the Integrating Matrix 

The integrating matrix is a means by which a continuous function may be integrated 
with the use of a finite-difference approach. This numerical method is based upon the 
assumption that the function f(x) may be represented by a polynomial of degree r as 
given in the following equation: 

f(x) = aQ + a^x + + . . . + a r x r (1) 

Let fj denote the value of the function f(x) at the station x = x^ where 
i = 0, 1, 2, . . ., r. If the stations are assumed to be equally spaced, then x^ = Xq + ih 
where h is the spacing interval. Equation (1) may be expressed in the form of Newton's 
forward-difference interpolation formula (ref. 5) as 

f(x) = g(p) 

= f 0 + p- p Af 0 + p(p - l)A 2 f 0 + jr P(P “ 1)(P “ 2)A 3 f 0 

+ . . . + \ p(p - 1) . . . (p - r + l)A r f 0 (2) 

^ * 

where p is the dimensionless quantity 


The forward differences appearing in equation (2) are given by 

A k f 0 = (E - l) k f 0 (4) 

The shifting operator E is defined by E^ = f^ where j = 1, 2, . . ., r and i = 0 
for the interpolation formula in equation (2). 
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Since f(x) s g(p) and since from equation (3) dx = h dp, the integration of f(x) 
from xj to xj is given by 

f 3 f(x)dx = h O' g(p)dp ( 5 ) 

''xi •■'i 


In proceeding with the development of the integrating matrix, it is convenient to 
choose a specific polynomial. If the function is assumed to be approximated by a third- 
degree (r = 3) polynomial, then g(p) from equation (2) may be substituted into equa- 
tion (5) and integrated over each of the three intervals between xq and x 3 to give 


J 1 f ( x )d* = ^( 9f 0 + 19f l - 5f 2 + f 3) 
x 0 

(6a) 

f 2 f(x)dx = |-(-f Q + i 3fl + 13£ 2 - f 3 ) 
X 1 

(6b) 

y Xs f(x)dx = iL(f 0 . 5f j + 19f 2 + 9f 3 ) 

(6c) 


The function f(x) may be integrated over a large number of equal-length intervals 
by repeated use of equation (6b). Equation (6b) is chosen since it is applicable to the cen- 
ter interval and since it has the smallest associated error of the three equations. Thus, 
from equations (6), the numerical integration of f(x) over each of the n intervals from 
Xg to x n is given by 

J Xl f(x)dx = iL( 9 f 0 + 19f 1 -5f 2+ f 3 ) 

X 


o 

JJ 2 f(*)dx = ^(-f 0 + 13fi + 13f 2 - f 3 ) 


J 3 f(x)dx = + 1 3f 2 + 13f 3 - f 4 ) 


( 7 ) 


1 n_1 = ^(- f n-3 + 13f n-2 + 13f n-l " f n) 

x n-2 

£ n f(x)dx = |L(f n _ 3 - 5f n _ 2 + 19f n _i + 9f n ) 


n-1 
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The matrix jA^j is an (n + 1) square matrix. As noted earlier, this matrix is 
based upon a third-degree polynomial. In a similar manner, matrices such as this may 
be obtained by using polynomials of other degrees. In general, such matrices are identi- 
fied as [Ar] where r denotes the degree of the polynomial upon which it is based. 

The integral of f(x) from x Q to x t may be obtained by merely summing the 
integrals given for each interval from xg to x^. That is, 


f 1 f(x)dx = P f(x)dx + f 2 f(x)dx + . . . + f 1 f(x)dx 

J Xg J Xg J Xj X, , 

This equation is given in matrix notation for i = 0, 1, 2, 


x i 
1-1 

., n by 



( 9 ) 


(10a) 


(10b) 



(10c) 
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The matrix [^BJ is an (n + 1) square matrix having a lower triangular array of 
unit elements. 

Substituting equation (8a) into equation (10a) and replacing the matrix j^Agj by the 
general matrix [A^ gives 



The integrating matrix is now defined as 

[ I ]=[B][A r J (12) 

From this definition, equation (11) becomes 



which is the desired matrix expression. Thus, the premultiplication by the integrating 
matrix of a column matrix defining the function f(x) yields the numerical integration of 
f(x) from Xq to x^ where i = 0, 1, 2, . . ,,n. Note that the number of intervals n 
is equal to or greater than the degree r of the polynomial representation. 

As stated in the Introduction, integrating matrices based upon polynomials of seventh 
and lower degree have been determined. Appendix A gives the matrix |^A^j for 
r = 1, 2, . . .,7. 

In reference 6 an integrating matrix is developed by a different approach. The 
integrand is approximated over two intervals by passing a parabola through three equally 
spaced points. Expressions for the coefficients of the equation of the parabola are deter- 
mined by solving a set of three simultaneous equations. Then the integrating matrix is 
derived through a series of matrix operations. The integrating matrix obtained is the 
same as that given by the product |j3j Ja^J of equation (12) for the specific case of 
r = 2 and n = 10. In reference 6 the integrating matrix is used to evaluate the integrals 
of an "integral series" solution, which is similar to the matrizant method, instead of as 
an operator by which the differential equations may be directly integrated. 

By using Newton's forward-difference interpolation formula, the integrand may be 
represented conveniently by polynomials of any degree. The interpolation formula avoids 
the solving of a set of simultaneous equations since it implicitly contains the polynomial 
coefficients in terms of the values of the integrand at the equally spaced stations. Also, 
when the integrating matrix presented herein is applied, any number of stations may be 
chosen (so long as n ^ r) since the integrating matrix was developed in general for any 
value of n. 
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Differential Equations and Boundary Conditions 

The derivation of the differential equations of motion, along with the statement and 
discussion of the assumptions made, is given in appendix B. A segment of displaced 


z 



Sketch a 


blade is shown in sketch a. The equations of motion for a freely vibrating propeller blade 
having a rotational velocity of £2 are 



where the variation of the axial tensile force T with respect to the blade longitudinal 
coordinate x is given by 

— + f2 2 mx = 0 (16) 

dx 

Also, in equations (14) and (15) v and w are the in-plane (of rotation) and out-of-plane 
displacements, respectively; m is the distributed mass; w is the natural vibration 
frequency; E is the elastic modulus; I yy and I zz are the moments of inertia of area 
about cross-sectional body axes which are parallel and perpendicular, respectively, to the 
plane of rotation when the blade is in the undeformed position; and I yz is the product of 
inertia. These relations are applicable to twisted propeller blades having a nonuniform 
and unsymmetrical cross section. 
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The propeller blade is assumed to have cantilever boundary conditions. Also, the 
axial force at the tip of the propeller blade is zero. The boundary conditions are listed 
in equations (17) to (21), where the subscripts 0 and n indicate that the associated 
quantity is evaluated at the fixed end x = Xq or at the free end x = x n . 



T n = 0 (21) 

Matrix Integration of the Differential Equations 

Since the differential equations are applicable for all values of x, they may be 
written for each of a chosen set of equally spaced stations defined by Xj = x Q + ih where 
i = 0, 1, 2, . . ., n. Thus, each of equations (14) and (15) determines a set of (n + 1) 
equations. These two sets of equations may be combined into a single matrix equation by 
defining 
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( EI yy)o 



o 



<Fy*)o 




(Etyz)i 

■v 


0M„ 



(22d) 


\ 

0 

0 




1 



(22e) 


where all off-diagonal elements not shown are zero. Also, each partition of is a 

diagonal matrix. By use of the partitioned matrices defined in equations (22), equa- 
tions (14) and (15) may be expressed in matrix notation as 



where £l^ is the unit matrix and {C T}> is a null column vector. 

Since the coupled displacements of the propeller blade lead to a partitioned-matrix 
equation, the matrix operator necessary to integrate this matrix differential equation is 
given in terms of the previously developed integrating matrix by 
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(24a) 



(24b) 


Because the integration of the differential equations is to be performed in matrix 
notation, a clearer insight into the process may be gained by considering first the simple 
equation d£/dx = f(x,£). Integrating this equation from xq to xj gives 




+ C 1 


(25a) 


where the constant of integration Cj is equal to and may be determined from a 

boundary condition. The subscript i indicates that the function £ is evaluated at 
x = Xj where i = 0, 1, 2, . . n. In equation (25a) the integral, whose integrand is a 
function of the unknown dependent variable £, may be expressed by some numerical 
method in terms of the integrand's unknown values corresponding to x = Xj where 
j = 0, 1, 2, . . ., n. Thus, after the numerical integration is performed, equation (25a) 
may be rewritten as 


?i 



(25b) 


which gives as a function of the x-j and ^ values. Since i = 0, 1, 2, . . .,n, 
equation (25b) represents a set of (n + 1) simultaneous linear algebraic equations which 
may be solved for the values to give the numerical solution of the differential 
equation. 


The numerical solution just described is essentially that which is effected when the 
integrating matrix is applied to a first-order differential equation. However, when the 
integrating matrix is used, the set of (n + 1) equations is expressed in a matrix equation 
and all the numerical integrations are performed by a single matrix operation. The 
integrating-matrix approach also allows the numerical solution of higher order differen- 
tial equations to be developed in a compact and orderly fashion. For the propeller blade 
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problem a matrix equation, which is analogous to equation (25b) and which relates the 
displacements of stations along the propeller blade, is obtained. 


The second term of equation (23) may be integrated numerically by pre multiplying 
the term by the integrating matrix pf], Thus, the integration of the equation yields 

- MsW -HpC 1 !! + n2 C c ]1 C m ] O) + ( K x> - (?) < 26 > 


The constant-of-integration matrix is given in general by the 2(n + 1) column matrix: 


r k ^ 

k jl 


k. 


jl 


<W- V- 


]2 


k. 


]2 


^ k i 2 ^ 


(j = 1, 2, 3,4, 5) (27) 


Equation (16), which is the force equilibrium relation for the X-direction, may be 
written in matrix notation as 


where 




(28 a) 


(28b) 
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Also, the second term of equation (26) may be expressed as 

El s ® - [eM] «• 

Substituting equation (28a) into equation (29) and then substituting the result into equa- 
tion (26) leads to 

- H®) - 



(30) 


Integrating this matrix differential equation by again operating with the integrating 
matrix gives 




- h 2 [" 2 w + h o> * rac^o + ( k 2> = c°> < 3 » 

Equation (31) may be premultiplied by the inverse of the stiffness matrix [sj and 
then integrated twice. The result after each of the successive integrations is given in the 
following equations: 

jl { „> - H[s ]-‘H & - 

- + n2 [ c ]] h - racj'raiw 

+ P]W' 1 {K 2 } + {K3} =» (32) 
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{cp} - [i] 2 [s] _1 [t]{<p} - fi 2 [l] 2 [s] _1 [T] [m] [xj{<0> 

- + ^ 2 H] H © + P] 2 H _1 H{ k i} 

+ [l] 2 [s] _1 {K2} + P]{K 3 } + (K 4 ) = 0 


(33) 



The operator matrices jjDqj and jD n J , which are used in evaluating the constant - 
of-integration matrices, are defined as 





1 

0 0 

0 

1 

1 




1 

0 

0 

1 

All zero 



1 

1 

1 

1 

1 

1 

1 

elements 


_ 

1 

0 

0 

1 

1 



- 




1 1 

0 0 

0 



All zero 


1 

• 1 

0 

0 



elements 


1 I 
1 1 

1 

1 


- 



1 1 

0 

0 


0 

0 0 

1 

1 

1 




0 ■ 

- - o 

1 

1 

1 

All zero 



1 

1 

1 

1 

1 

elements 


] - 

0 ■ 

0 

1 

1 

1 



1 




; o 

0 0 

1 



All zero 


! o ■ 

0 

1 



elements 


i i 
i i 

1 

1 

1 

1 





1 0 • 

0 

1 


(34a) 


(34b) 


In these matrices, each partition is an (n + 1) square matrix. 

The constant-of-integration matrix is evaluated by premultiplying equa- 
tion (26) by p) n J applying the boundary conditions given by equations (20) and (21). 

Since JjD^j > the premultiplication leads to 

{ K l} = [ D n] + ^ 2 [C]] [m| {</>} (35) 

Substituting this expression for pK-Q, into equation (31), premultiplying the 
resulting equation by , and applying the boundary conditions of equations (19) and 

(21) yields 

{K 2 } = SJ 2 [D n ] [l] [m] [x] {V} + [D n ][i] 2 p[l] + S-i 2 [c]J[m] 

- WDGIpJDG^ 1 ] + n2 [ c ]] H < f } (36) 
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Pre multiplying equation (32) by , applying the boundary conditions given in equa- 

tions (18), and noting that the first row of the integrating matrix H consists only of 
zero elements gives 

(K 3 ) = (0> (37) 

Similarly, the premultiplication of equation (33) by jjDqj and the application of the 
boundary conditions of equations (17) gives 

(K 4 ) = {0} (38) 

A relation for the normal axial force T may be obtained from equation (16). 
Integrating this equation numerically gives 


where 


{t} + + (k 5 ) . { 0 } 


(39a) 


<?>- » n - 

Tn 


Premultiplying equation (39a) by jjD^j and applying the boundary condition of equa- 
tion (21) gives 

{K 5 > = -n 2 [D n ][I][m][ X ]{l} 

If the matrix QfJ is defined as 

then equation (39a) may be expressed by 

(T> = n 2 [F][m][x]{l} 
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Equation (33) requires that the matrix representing the normal axial force T be 
a diagonal matrix; however, the column matrix of equation (42) results since the inte- 
grating matrix is used to obtain T. To transfer the elements Tj of the column matrix 
to a diagonal matrix, let ^PJ be a diagonal matrix whose diagonal elements are the 
same as the elements of the corresponding rows of the column matrix Qf] [\qj [x] {Y} . 
Therefore, 

[ T ] = fi 2 [p] (43a) 

where 


Y>] = diag[F][m][x]{l} 


(43b) 


The substitution of equations (35), (36), (37), (38), and (43a) into equation (33) yields 


[0 + X*] 2 [ s ]’ 1 [-KI + BBH - BYM 

Dividing equation (44) by w 2 and rearranging gives 

0W = »[»]W 

where 

00 * 

H = [i] ♦ rfCtt*]- 1 [-W - B HH - B 2 B y 


(44) 

(45a) 

(45b) 

(45c) 

(45d) 


Equation (45a) is the desired eigenvalue problem which results from the matrix integra- 
tion of the differential equations. The solutions of the eigenvalue problem define the 
natural vibration frequencies and the associated modal functions. 

Appendix C gives an expression for the stiffness matrix [jQ in terms of the 
moments and product of inertia about a pair of general orthogonal cross-sectional axes. 
Also, by utilizing the unique properties of the stiffness matrix, a relation is developed in 
this appendix which expresses explicitly the inverse of jjf) and thus avoids the inver- 
sion by usual computational methods. 
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Solutions of the Eigenvalue Problem 
The eigenvalue problem defined by equations (45) may be rewritten as 


where 


[n] {V(s)} = Ms) {<Ms)} 

[N] = 


(46 a) 
(46b) 


and (s) is used to denote the solution corresponding to the sth mode of vibration. Note 
that the matrix [h] is nonsingular and that the matrix [nJ is unsymmetrical. 


A close examination of the product £hJ ^jjcTj shows that the matrix (jsfj is of 
rank 2n since all the elements of the first and (n + 2)th rows are zero. Since the order 
of [n] is (2n + 2), there are at least two zero eigenvalues (A = 0). This may be proved 


by the expansion of the characteristic equation [~nJ - a£i^| = 0. 
sponding to the zero eigenvalues have no physical significance. 


The modes corre- 


A discussion of the many various methods of solving eigenvalue problems such as 
that of equations (46) is beyond the intent of this paper. However, two methods which have 
been used to obtain the solutions of equations (46) are explained briefly in the following 
sections. 


Solution using the QR transformation.- This method of solution is presently 
employed in the most general eigenvalue subroutine used by the computer center at the 
Langley Research Center. First, the matrix [nJ of equations (46) is reduced to upper 
Hessenberg form by elementary similarity transformations (refs. 7 and 8). Next, the 
similarity transformations, known as QR transformations (refs. 7, 9, and 10), of Francis 
are used iteratively to reduce further the matrix to an upper triangular matrix whose 
diagonal elements are the eigenvalues. The eigenvectors corresponding to the real 
eigenvalues are computed by using the inverse iteration method of Wielandt (as discussed 
in ref. 10). Since [V] is unsymmetrical, there is the possibility of complex eigen- 
values; however, complex eigenvalues have not been encountered in the use of this solu- 
tion method which does determine any imaginary components. 

The eigenvectors of the check cases on this subroutine were accurate to 10 sig- 
nificant figures. Thus, it is assumed that the errors associated with the numerical 
examples of this paper arise from the numerical integration approximations rather than 
from the solution of the eigenvalue problem. 

Solution using swee ping and iteration.- The sweeping and iteration procedure is a 
well-known and often-used method of solving eigenvalue problems. Iteration is a process 
by which an eigenvector is obtained by repeated premultiplications of a trial vector by a 
coefficient matrix. Applying the iteration process to equations (46) gives the first mode 
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of vibration. The higher modes are obtained by using a sweeping matrix to remove the 
lower mode components from the trial vector. The sweeping matrix is based upon the 
orthogonality relationship, which is derived in appendix D, between the modes. The 
development of the sweeping matrix for this problem is given in appendix E. Also, a 
proof of the convergence of the iteration procedure is given. 

NUMERICAL EXAMPLES 

The results of two natural vibration problems are presented to give an indication 
of the accuracy of the analysis and to show the effects of the variation of certain param- 
eters. The numerical examples considered are a typical propeller blade and a nonro- 

\ 

tating, cantilevered, uniform beam having lateral displacements in orily one direction. 

To verify the applicability of the differential equations of motion to a practical problem, 
the natural vibration frequencies of the propeller blade were determined numerically and 
compared with experimental data. The uniform-beam problem is analyzed to determine 
how the accuracy is affected by the choice of the integrating matrix and by the number of 
stations, as well as to substantiate further the accuracy of the presented method of 
solution. 

In appendix F, the differential equation of motion for the nonrotating beam having a 
symmetric cross section is integrated and formulated into an eigenvalue problem. This 
eigenvalue problem, along with the one of the propeller blade, was programed for solution 
by digital computer. For the propeller blade, the inverse of the stiffness matrix as given 
by equations (C6) and (C13) in appendix C was programed. The solutions to the eigen- 
value problems were obtained by using the previously described computer subroutine that 
utilizes the QR transformation. 


Propeller Blade 

The propeller blade selected for analysis is the WADC S-5 scale model of refer- 
ence 11. This blade was chosen since this reference gives a structural description suf- 
ficient for the numerical solution as well as experimental data for the natural vibration 
frequencies. The propeller blade is in effect cantilevered at 6 inches (0.1524 m) from 
the center of rotation and the tip of the blade is at a radius of 24 inches (0.6096 m). The 
ratio of the blade thickness to chord length varies from a value of 0.064 at the cantilever 
radius to a value of 0.021 at the tip. The stiffness characteristics are expressed in 
terms of the moments of inertia about the principal axes whose orientation with respect 
to the plane of rotation is given by the twist angle j3 (see sketch f in appendix C). The 
variation of j3 from the root to the tip of the blade is about 40°. In the experimental 
program, tests were made for various angles of pitch. The pitch settings were defined 
by the values of /3 as measured as x = 0.7 5R where R is the radius from the center 
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of rotation to the tip of the blade. With the exception of the rotational speed, all the input 
data necessary for numerical solution are given in table 1 for the pitch setting of /3 = 0 ° 
at 0.75R. Of course, for a different setting, the only change in the input data is a change 
in the values of /3 by an amount equal to the setting. 

In order to compare numerical results with test data, solutions were computed for 
cases corresponding to the pitch settings and rotational speeds of the experimental inves- 
tigation. The first- and second-mode natural frequencies, but not mode shapes, were 
determined experimentally in an evacuated chamber for various rotational speeds at pitch 
settings of 0°, 20°, and 40°. The numerical solutions were obtained by using 10 stations, 
which correspond to nine 2-inch (0.0508-m) intervals, to describe the cross-sectional 
properties of the propeller blade. Also, a fifth-degree integrating matrix was used in 
computing the solutions. The terminology "fifth degree" denotes the integrating matrix 
based upon a fifth-degree polynomial (r = 5) . 

The experimentally and analytically determined natural vibration frequencies are 
given in table 2 for the various cases. This table also gives the percent error which is 
the difference between the experimental and computed values expressed as a percentage 
of the experimental value. These results are presented in the manner typical of 
propeller-blade analyses in figure 1, which shows graphically the natural frequencies as 
a function of the rotational speed. The dashed lines (labeled IP, 2P, and 3P) of the fig- 
ure relate the natural frequencies as dimensionless multiples of the rotational speed. 

The intersections of the frequency curves with these dashed lines define the critical 
rotational speeds (i.e., the speeds at which excitation is most likely to occur). It is 
observed that the existence of a crossing of the IP line is dependent upon the pitch of the 
propeller blade. It is also noted that the pitch has very little effect on vibration fre- 
quencies of modes above the fundamental. 

A comparison of the modal displacements is not possible since the mode shapes 
were not determined in the experimental investigation. However, for the sake of pre- 
senting examples of typical propeller modes and showing how the mode shapes are affected 
by rotational speed, figure 2 illustrates the three lowest computed modal functions at 
rotational speeds of 0 and 6016 rpm with the blade having a pitch angle of /3 = 20° 
at 0.75R. 


Uniform Beam 

Numerical solutions for the natural vibration characteristics of a cantilevered uni- 
form beam were computed by using various integrating matrices and various numbers of 
stations. In order to simplify the comparison of the numerically determined natural fre- 
quencies with the frequencies of the exact solutions as given by reference 12, the inputs 
m and h were chosen to be equal to unity and El was set equal to n^. Thus, the 
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quantity El/mZ^ j s always equal to unity since the length l is equal to nh. For such 
a beam, the exact natural vibration frequencies for the first eight modes are as follows: 


Mode 

Vibration frequency, o>, 
rad/sec 

1 

3.5160154 

2 

22.034492 

3 

61.697214 

4 

120.90192 

5 

199.85953 

6 

298.55553 

7 

416.99079 

8 

555.16525 


Figures 3 to 9 present the percent errors, as defined by * — computed ^ — exact) — 

exact 

of the natural frequencies computed with the use of integrating matrices of degrees one 
to seven, respectively. Each figure gives the error percentage of the frequencies for tl 
first eight modes as a function of the number of stations used in the computations. The 
dashed curves in the figures denote percent errors having a negative sign. 

A small portion of the generated modal data is presented in order to give an indi- 
cation of the accuracy of the various computed mode shapes and to show how the accura 
is affected by the degree of the integrating matrix and by the number of intervals. 

Tables 3 to 6 give a comparison of the modal displacements computed by using the vari 
ous integrating matrices with those of the exact solution. For the computation of the di 
placements, both 25 intervals (n = 25) and 10 intervals (n = 10) were used, and the dis- 
placement at the free end of the beam was taken as two units. However, the data are 
presented only at x/Z increments of 0.2, which is adequate for comparison purposes 
and for indicating the accuracy of the computed modal functions. 

DISCUSSION 
Discussion of Accuracy 

The results presented in the previous section show that the governing differential 
equations are indeed applicable to propeller blades. Table 2 gives a comparison of com- 
puted natural frequencies with experimentally determined frequencies. The fact that the 
errors remain relatively small over the wide range of rotational speed and pitch setting 
indicates that the differential equations adequately describe the motion of the typical 
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propeller blade. The error percentage averages less than 2 percent with a maximum 
value of 3.6 percent. All the computed natural frequencies are greater than the experi- 
mental frequencies except for one case. This indicates that slight inaccuracies may 
possibly exist in the input data for stiffness and/or distributed mass. In addition to 
establishing the applicability of the governing equations, these results indicate that the 
presented integration method is an accurate numerical method for solving differential 
equations. 

An examination of figures 3 to 9, which give the percent error of the computed 
natural frequencies for a uniform beam, reveals the significance of the degree of the 
polynomial representation used in developing the integrating matrix. For a given number 
of intervals, it is seen that the use of the first-degree integrating matrix, which is based 
upon trapezoidal integration, results in errors much greater than those given by matrices 
based upon second- and higher degree polynomials. The errors of the frequencies com- 
puted with the second- and third-degree integrating matrices are of the same magnitude, 
with the third-degree matrix being slightly more accurate. The fourth-degree integrating 
matrix yields accuracy which is significantly better than that of the third-degree matrix. 
For an error limit in the frequencies of about 1 percent, which is satisfactory for most 
engineering applications, the figures show that the fourth-, fifth-, sixth-, and seventh- 
degree integrating matrices each require about the same number of intervals. When 
these integrating matrices are used, it is seen that an accuracy of 1 percent will result 
if the number of intervals is slightly greater than twice the number of the highest mode 
to be determined. However, to obtain very small errors in the frequencies, the number 
of intervals required appears to decrease as the degree of the integrating matrices 
increases. 

The seventh-degree integrating matrix (fig. 9) yields errors of 0.001 percent when 
the number of intervals is approximately four times the number of the mode. The results 
for this example indicate that the solution method would be very useful in applications 
where extreme accuracy is desired. For example, the first-mode frequency had an error 
of only 0.000014 percent when the seventh-degree integrating matrix with seven intervals 
was used. 

The curves of these figures show that the solution is apparently stable as the num- 
ber of intervals increases. There is an oscillation (change in sign) in the value of the 
percent error for a few of the curves (see figs. 5 to 9); however, in all cases the numeri- 
cal solution appears to converge, as expected, to the exact solution as the number of 
intervals increases. 

Table 3 shows that all the computed first-mode displacements are very accurate. 

In fact, a maximum error of 0.00001 in the modal displacements is given by all the inte- 
grating matrices except the first-degree matrix when 25 intervals are used and the 
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first- and second-degree matrices when 10 intervals are used. The higher modes 
(tables 4 to 6) show that the accuracy increases with the degree of the integrating matrix. 
When the seventh -degree matrix and 25 intervals are used, the maximum error in the com- 
puted modal displacements continues to be 0.00001 for the higher modes. As expected, 
the accuracy is significantly reduced for the higher modes that are computed with 10 inter- 
vals. However, with only 10 intervals the seventh-degree integrating matrix yields results 
for the fourth mode which are satisfactory for most engineering problems. The computed 
fourth-mode characteristics have a maximum error in the modal displacements of less 
than 0.01 and an error in the natural vibration frequency of 0.52 percent. 

The uniform-beam vibration problem demonstrates that in general, the accuracy of 
the computed vibration characteristics increases with the degree of the integrating matrix 
and with the number of stations. Also, it shows that for a given accuracy and a given 
mode the number of stations required is small relative to the number needed for other 
numerical methods. Obviously, the higher degree integrating matrices are to be preferred 
since the same accuracy can be obtained with a fewer number of stations, which means 
less computational time. Integrating matrices based upon polynomials of degree greater 
than seven may be developed, but the presented results indicate that the improvement in 
the accuracy becomes less with each increase in the degree of the polynomial. 

The accuracy of the uniform-beam solutions is thought to be typical of the accuracy 
which may be expected for the propeller-blade problem and for other beam vibration 
problems which use the integrating matrix to obtain solutions. The uniformity of the 
example beam is not thought to be a factor that enhances significantly the accuracy of the 
method since it is the modal function which is first integrated numerically. 

The integrating matrix yields exact integrals when the integrand is of degree equal 
to or less than the degree of the polynomial upon which the integrating matrix is based. 
However, when the degree of the integrand is greater than the degree of the polynomial, 
the polynomial representation approximates the integrand, and thus, the integrals given 
by expressions such as equations (6) become numerical approximations. The error asso- 
ciated with each integral approximation may be determined by expressing the integral as 
a Taylor's series and then substituting finite -difference approximations for the deriva- 
tives. It can be shown that the errors of the approximations used in developing the inte- 
grating matrices are o(h r+2 ) where r is the degree of the assumed polynomial. For 
example, the error of each of equations (6) is of order h^. 

The QR transformation method, which obtains all the eigenvalues at once and does 
not depend on the eigenvectors, is much preferred over the sweeping and iteration tech- 
nique. The iteration process yields accurate eigenvalues, but the slight errors which 
occur in the eigenvectors usually prevent the determination of solutions above the third 
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or fourth mode. A possible reason for this behavior is the fact that the orthogonality 
relationship is derived from the differential equations rather than from the matrix 
equation. 


Other Applications of Solution Method 

The integrating matrix may be used to solve various other boundary-value problems 
in addition to determining the lateral vibration characteristics of propeller blades and 
beams. For example, the author has used the integrating matrix to develop solutions for 
the longitudinal vibration characteristics of beams, the vibration characteristics of beams 
having coupled lateral and torsional motions, the buckling characteristics of columns, 
and the lateral displacements of beams under applied loading. When the integrating 
matrix is used to solve for the natural vibration characteristics of beams, an eigenvalue 
problem results in which the dominant eigenvalues correspond to the lower frequencies 
without the need for determining the inverse of the coefficient matrix. Also, for the 
nonhomogeneous-beam problem, the solution for the displacements is given directly with- 
out having to solve a set of simultaneous equations. For coupled differential equations 
which are of different order, such as those of a beam having bending and torsional dis- 
placements, the solution may be formulated by integrating each equation separately and 
;hen combining the resulting matrix equations into a single partitioned-matrix equation. 

In addition to boundary-value problems, the integrating matrix may be applied to 
nitial-value problems with the solution being developed in the same manner. Although 
;he usefulness of the integrating-matrix method for solving typical engineering initial - 
/alue problems has not been demonstrated, simple examples, whose exact solutions are 
mown, show that the method yields results which are extremely accurate relative to those 
}f other methods (e.g., the Euler and center-slope methods) for the same spacing interval. 

As an example, the solution to the first-order equation + w = 0 with the initial con- 
dition of w(x=0) = Wq is found very simply to be 

{*> = *o[i 1 ] + H] 0> 

where {w} is of the same form as <^w} which is defined by equation (F3b). After the 
inverse is evaluated, each w^ is related to Wq by a constant. To obtain the solution 
w. for i > n, w n becomes the initial condition and replaces w Q in the solution. This 
process may be repeated indefinitely to obtain WqjWj, . . ., w n , . . ., w 2n > • • •> w 3n , 
and so forth. It is noted that the inverse need only be computed once. Also, the solution 
may be obtained by calculating only the values w n ,w 2n ,w 3n , . . ., if the intermediate 
values are not desired. For this case, the spacing interval in effect becomes nh. By 
choosing the integrating matrix corresponding to r = n = 2 with h = 0.1 and letting 
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Wq = 1.0, the solution given in the following table is computed and compared with the 
exact solution: 


X 

w 

Exact 

Computed 

0 

1.0 

1.0 

.2 

.818731 

.818731 

.4 

.670320 

.670321 

.6 

.548812 

.548812 

.8 

.449329 

.449330 

1.0 

.367879 

.367880 


Thus, it is seen that a very accurate solution results from using only a second-degree 
integrating matrix. Upon first inverting the 3x3 matrix of the solution formulation, each 
Wj was obtained simply by one multiplication. 

The integrating matrix may also be applied to initial-value problems which are non- 
homogeneous and higher order. Furthermore, it is not necessary that the spacing interval 
be constant as was the case in the preceding example. For the nonhomogeneous equation 

■T— — + w = f(x), the solution at each station is found to be given by 


{w} = ^[l] +[l]J j^ 0 {l} + 


where Wq may be replaced as before by w n , W2 n , and so forth. Examination of the 
last row of this matrix equation shows that the solution at every nth station is given by 


n 

*(j+l)n = a ^jn + 2 b i f jn+i 

i=0 

where j = 0, 1, 2, . . ., and a and bj are constants resulting from the matrix 
multiplications. 


Discussion of Method 

The integrating-matrix method of solving differential equations is a straightforward 
approach (e.g., see appendix F). That is, the differential equations are integrated very 
simply by using the integrating matrix as an operator, and then the constants of integra- 
tion are evaluated from the boundary conditions. In the presented propeller-blade prob- 
lem, the development of the solution is complicated considerably by having coupled 
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differential equations and by the centrifugal -force effects. For boundary- value problems, 
such as those mentioned in the previous section, the solution is formulated very quickly 
by this method. The solutions are expressed concisely and are readily programed by 
using standard subroutines in most instances. In obtaining the numerical solutions, the 
inputs to the computational programs are given simply by the local values of the cross- 
sectional properties; the generation of stiffness or influence coefficients is avoided. The 
computed solutions are accurate because high-degree polynomials are used to approxi- 
mate the functions which are integrated. 

As noted in the Introduction, the method is applicable to continuous problems rather 
than to discrete or lumped mass systems. Since the polynomial representations are of 
finite degree, the method is best applied to problems in which all functions are smoothly 
varying. However, with discretion the method may be applied to, for example, the prob- 
lem of a beam which has an abrupt change in its cross-sectional properties such as the 
stiffness or the distributed mass. Problems having abrupt changes require a smaller 
interval size. For such problems it is prudent to obtain a second solution by decreasing 
the interval size and/or by increasing the degree of the integrating matrix. If either of 
these changes alters the solution appreciably, then it may be concluded that the interval 
size was too large in the first solution. The varying of the interval size without changing 
the solution indicates that the solution is valid. 

CONCLUDING REMARKS 

A numerical method for determining the natural vibration characteristics of rotating 
beams, such as propeller blades, is presented. By applying the integrating matrix as an 
operator, the derived differential equations are integrated and result in a matrix eigen- 
value problem. Two methods of solving the eigenvalue problem are described, one of 
which utilizes the QR transformation and inverse iteration. The other method is a 
sweeping and iteration procedure that requires a special sweeping matrix based upon the 
derived orthogonality relationship. The stiffness properties may be transformed to any 
pair of orthogonal cross-sectional axes, and an explicit expression for the inverse of the 
stiffness matrix is derived. 

The computed natural frequencies for atypical propeller blade were compared with 
measured frequencies. In addition, an error analysis of the numerical method was made 
by comparing the computed vibration characteristics of a uniform beam with exact 
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solutions. These two numerical examples showed that the method of solution yields very 
accurate results. A discussion of the possible application of the solution method to other 
problems is included. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., September 18, 1970. 
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APPENDIX A 


INTEGRATING MATRICES 


The integrating matrix ^ij was defined by equation (12) as 

H-HM 

where the subscript r denotes the degree of the assumed polynomial upon which the 
integrating matrix is based. The matrices are as follows for r = 1, 2, . . .,7: 
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APPENDIX A - Concluded 
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APPENDIX B 

DERIVATION OF THE DIFFERENTIAL EQUATIONS OF MOTION 

In this derivation of the equations of motion for a freely vibrating propeller blade, 
the propeller blade is idealized as a rotating beam having a nonuniform and unsymmet- 
rical cross section. This beam is depicted in sketch b in the unstrained position. The 
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Sketch b 


centroidal axis, which is assumed to be straight, lies along the X-axis of the rotating 
XYZ orthogonal coordinate system, which has its origin fixed and is rotating about the 
Z-axis with an angular velocity of S2. The elastic axis of the beam is assumed to be 
coincident with the centroidal axis; in other words, the bending and torsional motions are 
assumed to be uncoupled. This is a sound assumption for the lower vibration frequencies 
of propeller blades because the proximity of the elastic and centroidal axes renders small 
coupling. The comparison of experimental vibration data of typical propeller blades with 
the numerical results of this and other analyses shows that the errors induced by 
neglecting torsional motion are very small for the lower bending modes. Also, the 
secondary effects of shear deformation and rotary inertia are considered negligible for 
propeller blades since the cross-sectional dimensions are small compared with the length. 
In addition, the cross-sectional properties of the beam are assumed to be continuous func- 
tions with no abrupt changes, and the lateral displacements are assumed to be small. 

The inertial forces acting on a differential segment of the rotating beam are deter- 
mined from the acceleration vector of the centroid of the segment. The effect of the lon- 
gitudinal displacement of the centroid on the inertial forces is small and is assumed 
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Sketch c 

negligible. The position vector (see sketch c) to the centroid of a general differential 
segment of the strained beam is given by 

r = xi +vj + wk (Bl) 


where x is the body-axis coordinate in the X-direction for the unstrained beam. Also, 
v and w are functions of x and represent the elastic displacements of the centroidal 
axis in the Y- and Z-directions, respectively. The absolute angular velocity of the 
rotating system is assumed to be constant and is given by u> = fik . Differentiating equa- 
tion (Bl) twice with respect to time gives the acceleration vector as 


a = 



(B2) 


It can be shown that the Coriolis acceleration term may be neglected. The natural 
frequencies which are of particular interest to propeller -blade analyses have the same 
order of magnitude as the rotational speed f2. If harmonic motion with the frequency 
being the same order as the rotational speed is assumed, the Coriolis acceleration term 
is obviously negligible when compared with D 2 x since the amplitude of v is much less 
than the magnitude of x for small displacements. 


Since the acceleration components are known, the reversed effective inertial forces 
acting in the X-, Y-, and Z-directions on a differential segment of length dx are 

fi 2 xm dx, (f2^v - -2-^m dx, and - m dx, respectively, where m is the distributed 

V at 2 / at 2 


mass of the beam. 
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Sketch d 

The inertial forces, internal forces, and internal moments lying the the XZ-plane 
are shown acting on a differential beam segment in sketch d where T is axial tensile 
force, V is shear force, and M is bending moment. The summation of forces in the 
X-direction and in the Z-direction gives the two following equations: 


— + fi 2 mx = 0 (B3) 

ax 


3V Z 

dX 


m ^-=0 


at 2 


(B4) 


After higher order terms are dropped, the summation of moments in the Y-direction 
yields 


V 


Z = 


9M y 

+ 

Ox 



(B5) 


Substituting equation (B5) into equation (B4) gives 


9 2 M„ 





m 



= 0 


(B6) 
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Sketch e shows the forces and moments acting on the segment in the XY-plane. By 
proceeding in the same manner as before, the summation of forces in the Y-direction and 
the summation of moments in the Z-direction yields 


3V Y 
9x 


+ 



m = 0 


(B7) 


V Y =- 


9M Z + T 9v 

ax ax 


Upon substitution of the expression for Vy, equation (B7) becomes 




m = 


0 


(B8) 


(B9) 


Relations for the bending moments My and Mg in terms of the displacements 
v and w may be determined by integrating over the cross-sectional area the products 
of the longitudinal stress and each of the body-axis coordinates of which the stress is a 
function. It is necessary to determine first an expression for the longitudinal strain over 
the cross section. The cross-sectional body-axis coordinates are y and z which are 
parallel to the Y- and Z-axes, respectively, when the beam is in the unstrained condition. 
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Assuming that plane sections remain plane and making small-angle assumptions 
for dv/dx and aw/ ax, which define the slope of the centroidal axis, gives the longi- 
tudinal displacement of a point defined by y and z within a cross section as 


u = u 




(BIO) 


where u is the elastic displacement of the centroidal axis in the X-direction. Since the 

longitudinal strain is defined by e = — , equation (BIO) leads to 

9x 


e = 


_9u _ 8 2 V 

8x 3x2 



(Bll) 


From Hooke's law, the longitudinal stress is merely the product of the strain and thf 
modulus of elasticity since the normal stresses are zero. If the elastic modulus E is 
assumed to be constant over the cross section, the bending moments are given by 

My = E \ ez dA 
J A 


M Z = 



(B12) 


Substituting equation (Bll) into equations (B12) and integrating over the cross-sectional 
area yields 


M v = -El 


yy 


a 2 w 

9x2 


- EI 


y Z 


9 2 v 

ax 2 


(B13) 


M Z = EIyz-^f +EI ZZ ^| 

ax 2 ax 2 

where the moments and product of inertia are defined by 


‘yy = 5 

* 

z 2 dA 

A 

‘zz = j 

y 2 dA \ 
A ' 


\ 

yz dA 

A J 


(B14) 


(B15) 
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The substitution of equations (B13) and (B14) into equations (B6) and (B9), respec- 



An examination of equations (B16) shows that the two displacements v and w 
are coupled by the product of inertia I yz of the unsymmetrical cross section. There- 
fore, for I yz * 0, any deflection in the Y-direction has associated with it a deflection in 
the Z-direction. The static coupling between the lateral displacements may be verified 
by using the bending-moment relations of equations (B13) and (B14) or by using directly 
the strain expression of equation (Bll). 

By assuming harmonic motion such that v = ve iwt and w = we iwt , the equations 
of motion become 


d 2 




d 2 w 

dx 2 



u> mw = 0 


(B17) 


d 


2 


dx^ 





0 


(B18) 


where u> is the natural vibration frequency. These two equations in conjunction with 
the auxiliary relation of equation (B3), which may be rewritten as 

— + fi2 mx _ o (B19) 

dx 

are the governing differential equations of motion. 
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APPENDIX C 


TRANSFORMATION AND INVERSE OF THE STIFFNESS MATRIX 

It may often be desirable to express the stiffness matrix [jsj in terms of the 
moments and product of inertia about a pair of orthogonal cross-sectional axes other than 
the axes associated with the y,z body coordinates. As illustrated in sketch f, y' 

2 * 2 

y’ 
y 

Sketch f 



and z' are defined to be the coordinates of a reference system which is obtained by the 
rotation of the yz reference through an angle /3. The moments and product of inertia 
in the primed system are defined by 



(Cl) 


Also, the coordinate transformation relations are 


y* = y cos /3 + z sin /3 
z' = -y sin j3 + z cos /3 


(C2) 


By substituting equations (C2) into equations (Cl), noting that the Jacobian relating the 
differential areas in the primed and unprimed systems is equal to unity, and using the 
definitions of equations (B15), the moments and products of inertia of the two systems 
are found to be related by 
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2 .2 

IytyT = Iyy cos + I zz sin j3 " 2Iy Z SU 1 COS j3 

2 2 \ 

I z * z t - Iyy sin /3 + I zz cos /3 + 2Iy Z sin j8 cos /3 / 

Iy T z T = (iyy ” ^zz]®^ cos /3 + Iy Z (cos j3 ~ sin 
From equations (C3), the following matrix identity may be written: 

O'] ■ M T MM 

where 



M 



(cos /3)q 

(sin p) Q 


(cos /3) 1 

X 

X 

(sin 0). 

X X 

X 

N 


X 

(cos /3) n 

X 

(sin j3) n 

- 

-(sin 0) Q 

(cos 0) o 


-(sin p) l 

X 

X 

(cos 

V 


X 

X 

-(sin 0) n 

V 

V 

(cos j8) n 


(C3) 


(C4a) 


(C4b) 


(C4c) 
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and Q=f] is defined by equation (2 2d). Since the transformation matrix B is orthog- 
onal, its inverse is equal to its transpose; thus, from equation (C4a) the stiffness matrix 
in the yz reference may be expressed in terms of the stiffness values in the y'z' 
reference by 

[S]=[L][S'][L] T <C5, 

Also, taking the inverse of both sides of equation (C5) gives the transformation of the 
inverse of the stiffness matrix as 

<«> 

Since it is desirable to avoid the inversion of a matrix whenever possible, an explicit 
expression will be derived for the inverse of the stiffness matrix which appears in the 
eigenvalue problem of equations (45). Consider the matrix £a] given in equation (C7) 
which is of the same form as the matrices [V] and [s'] defined, respectively, by equa- 
tions (22d) and (C4b). However, in the development which follows it is not necessary for 
B to be symmetric as are []s[] and [s[] . 


V 


l 2>2„ 


1,1 




l n+l,l 


n+2,2 


a 

n+i,i 


a ' 
2n,n 


l l,n+l 


l 2 ,n+2 

7 s. 


i,n+i 


n,2n 


m+l ,n+l 


*n+2 ,n+2 


a n+i,n+i 


2n,2n 


(C7) 


In the above matrix [a] all off-diagonal elements in each of the four n-square partition 
matrices are zero and i = 1, 2, . . ., n. 

Notice that the cofactor of each zero element of [a] is zero since the minor of 
each such element contains two rows, as well as two columns, which are proportional. 

Let a., represent the cofactor of the element a. ., where j,k=l,2, . . ., 2n. Then 
the inverse of the matrix lAl is given by 
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M’ 1 “ TR 


H 


a 


1,1 


a 


2,2 


O'- . 


a 


a 


l,n+l 


a 2,n+2 


a- _ , . 
i,n+i 

s 


a; 


n,2n 


a 


n+1,1 


n+2,2 


Q! 


n+i,i 


a 


2n,n 


“n+l ,n+l 


a 


n+2,n+2 


a 


n+i,n+i 


a 


2n,2n 


(C8) 


The cofactor ^ in the upper-left partition of * is equal to the minor of 

the element a^ ^ with a sign. Since the minor of a^j is equal to the sum of the pro- 
ducts of the elements of the (n + i - l)th row (or column) of the minor and their respec- 
tive cofactors, the cofactor a. , may be expressed by 


a i,i ~ M i,i a n+i,n+i( -1 ) ^ ^ D i a n+i,n+i D i 


(C9a) 


where Mi i is the minor of the element a^ and Di is the determinant of the 
matrix [a] with the ith and (n + i)th rows and the ith and (n + i)th columns removed. 
In a similar manner, the cofactors of the other partitions are found to be 


“i,n + l = (- 1 )” +2lM l,n + i = = - a n«,i D ‘ 


n+2i-l. 


(C9b) 


% + i,i = <-« n+2lM n + i,i = = -\n + i D l 


n+2i-l . 


(C9c) 


a n+i,n+i = Mn +i , n+ i - 


2i, 


(C9d) 


The determinant of £a] may be expressed as the sum of the products of the ele- 
ments of the ith row and their respective cofactors. Thus, 
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. A J| + a i,n+i Q! i,n+i ( a i,i a n+i,n+i " a i,n+i a n+i,i) D i 


(CIO) 


From equations (C 9 ) and (CIO), it is seen that the elements in of equa- 

tion (C 8 ) are given by 


a. . 

1,1 _ 


~ a n+i,n+i 


M 


a. .a . . - a- -a 

1,1 n+i,n+i i,n+i n+1,1 


(Cl la) 


a i,n+i "^l+i.i 


|M 


a i,i a n+i,n+i “ a i,n+i a n+i,i 


(Cl lb) 


a 


-a. 


i,n+i 


n+i,i 

All a i,i a n+i,n+i - a i,n+i a n+i,i 


(Cllc) 


a 


n+i^n+i 


a, . 


a] 


a i,i a n+i,n+i ” a i,n+i a n+i,i 


(Cl Id) 


Substitution of the corresponding elements of equation (C 4 b) into equations (Cll) and 
then substitution of the results into equation (C 8 ) gives the inverse of [jS'J as 


M- 1 - 


EI z’z’ 


iEly, y ,El z , z , - El^ zt/o 


EI z’z’ 


\EI y , y ,EI z , z , EIy, z , 


n 


- E Vz> 


i E Vy EI z*z- - ei ;-z j Q 


• EI y , z- 


^ E 1 y f y fEI z TzT " EI y tz V n 


- E Vz' 


^Vy- E Vz : - e w 0 


~ EI y'z’ 


\ EI y'y’ EI z’z’ “ Ei ;- Z - 


n 


EI y’y’ 


L EI y»y- EI Z ' Z ^- E Vzy 0 


EI y’y' 


lEI y'y ,EI z’z' “ EI y’z’ /n 


(Cl 2) 


It is obvious that a similar expression may be written for QTj \ 
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If the angle 



/3 represents the orientation of a principal axis of inertia, then since 
reduces to the following diagonal matrix: 
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APPENDIX D 


DERIVATION OF THE ORTHOGONALITY RELATION 


A knowledge of the orthogonality relation between the natural modes is important 
to the analysis of many vibration problems. The orthogonality relation is often used to 
obtain the solutions of eigenvalue problems by iteration methods. Also, the relation is a 
very useful tool in modal series solutions of response problems. Since the displace- 
ments w and v are coupled, the conventional orthogonality relation for beams in 
bending is not valid. 

The differential equations of motion for the rth mode are 


dx^ 


El 


yy 


d 2 w(r) 

dx 2 


+ EIy Z 

* dx^ 


£[t - ct> 2 (r)mw(r) = 0 (Dl) 


- ^|t - [co 2 (r) + n^jmv(r) = 0 (D2) 

Multiplying equation (Dl) by w(s)dx and equation (D2) by v(s)dx, where w(s) and 
v(s) are the displacements of the sth mode, and then integrating both equations from Xq 
to x n gives 


dx2 


ET d 2 w(r) FT d 2 v(r) 



d 2 v(r) 

dx 2 


dx 


and 




mw(r)w(s)dx = 0 



d 2 w(r) 

dx 2 


+ EI ZZ 


d 2 v(r) 

dx 2 


dx 


(D3) 



mv(r)v(s)dx = 0 


(D4) 
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Performing the integrations in equations (D3) and (D4) by parts and applying the boundary 
conditions yields 



yy 


+ ei 

dx 2 


yz 


d 2 v(r) 

d 2 w(s) 

dx 2 

dx 2 


and 



dx - o> 2 (r) f n mw(r)w(s)dx = 0 

CLX CtX 

x o 


(D5) 


rt 


FT d 2 w(r) FT d 2 v(r) 

yz “S2 + EIzz '5T 


d 2 v(s) 


dx 2 


dx 


+ ( * n T Szill &SL dx - [o> 2 (r) + S2 2 ] f X " mv(r)v(s)dx = 0 
J v dx dx u 


(D6) 


In like manner, two additional relations similar to equations (D5) and (D6) may be 
obtained. First the differential equations are written for the sth mode as 


df_ 

dx 2 


FI d 2 w(s) FI d 2 v(s) 

% dx2 +EIyz lP“ 


d_ 

dx 


dw(s) 

dx 2 


F 


(s)mw(s) = 0 


(D7) 


dx z 


EI yz + El aWl . f [ T SZfcl] - [o, 2 (s) + o 2 ]mv(s) = 0 (D 8 ) 

_ dx 2 dx 2 J J 


Equations (D7) and (D8) may be multiplied by w(r)dx and v(r)dx, respectively, and then 
integrated from Xq to x n . Making the multiplications, performing the integration by 
parts, and applying the boundary conditions yields 


r 


VT ^ w(s) d v(s) 

Iyy dx2 + Iyz dx 2 


d^w(r) 


cbF 


dx 


J X " T ^ dx - u> 2 (s) J*” mw(r)w(s)dx = 0 


(D9) 
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Concluded 
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APPENDIX E 


SWEEPING AND ITERATION PROCEDURE 

The orthogonality relation derived in appendix D may by use of the integrating 
matrix be expressed in matrix notation as 


L Q J[l][m] [<?( r)] (<p(s)} = 0 


(r * s) (Ela) 


where the row matrix (_Qj, which consists of two partition matrices with (n + 1) elements 
in each, is given by 


Q = 


0 ... 0 1 < 0 


...01 


(Elb) 


An arbitrary vector may be expressed as a linear combination of the inde- 

pendent eigenvectors by 


(j&y - + c 2 + • • • + c 2n+2 + 2)} 


(E2) 


In equation (E2) the normalized eigenvectors jj)> correspond to the eigenvalues \(j) 
which decrease in magnitude as j varies from 1 to (2n + 2). An expression for the 
scaling constants Cj may be obtained by premultiplying equation (E2) by 

L Q J[l][m] [<p ( j)] and applying the orthogonality relation given by equation (El). Hence, 


c^ = 


LqJHH fed £0 


(E3) 


If the first (s - 1) modes are assumed to be known, the sth mode is obtained by first 
rewriting equation (E2) as 


{$} - c 1 - (cp( 2)} c 2 - ... - (cp ( s - 1)} c g _ 1 


= C g ^\p(s)y + c g+1 «J>(s + 1)} + . . . + c 2n+2 {cp( 2n + 2)} (E4) 

Upon substituting equation (E3) for c^, c 2 , . . ., c g _p equation (E4) becomes 

0<s)] {*} = c s {^(s)} + c g+1 (cp ( s + 1)} + . . . + c 2n+2 {y(2n + 2)} (E5) 
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where 


r*(s)i m &w] fr^LoJPlHM 

L q J P]H M < G p(i) ) L Q JP]H ^ (2) ) 


{y(s - 1 ))LoJP][ m ][^(s - 

LQj[l][m] [,(s - l)J{^(s - l)} 


(E6a) 


The square matrix (jKsj] is the desired "sweeping" matrix. For computational pur- 
poses, it is expressed conveniently by 


[*<s>] = [*( S - d] . <y<° - ^LqJCOHIj^ - »>] 

LqJP]W[^ s - nk** - d) 


(E6b) 


where is the identity matrix. 

Premultiplying equation (E5) by [Vj and substituting equation (46a) into the terms 
on the right-hand side of the resulting equation yields 


[n][*(s)J{$} = c g X(s) .(^(s)} + c g+1 X(s + 1) {^(s + lj) 
+ . . . + c 2n+2^^ 2n + + 


(E7) 


Defining [n(s)J as 

[N(s)] = [n][*(s)] (E8) 

and premultiplying equation (E7) by [n(s)J gives 

[n(s)] 2 {¥} = CgX(s) [n] |>(s3 ^(s)} + c g+1 X(s + 1)[N] [*(*)] {V(s + 1)} 

+ . . . + c 2n+2 X(2n + 2 )[n] QKs)] fy(2n + 2)} (E9) 


When the orthogonality relationship of equation (El) is applied and equation (46a) is sub- 
stituted into the resulting terms on the right-hand side, equation (E9) becomes 
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[n(s)] 2 {#} = c s X 2 (s){V(s)} + c g+1 X 2 (s + l){V(s + 1)} 

+ . . . + c 2n+2 X 2 (2n + 2){^(2n + 2)} (E10) 

The repeated premultiplication of equation (E7) by [n ( s)J gives 

[N(s)] k {¥) = c g X k (s) {<^(s)} + c g+1 X k (s + l){V(s + 1)} 

+ . . . + c 2n+2 X k (2n + 2) (cp(2n + 2)} (Ell) 

For k sufficiently large, equation (Ell) reduces to 

[N(s)] k {$} * c s X k (s) {cp( s)} (El 2) 


since X(s) > X(s + 1) > . . . > X(2n + 2). It is seen that the product obtained by repeated 
premultiplication of an arbitrary vector by jjN(s)J converges, as k increases, on the 

normalized eigenvector s)} times a scalar. Hence, the normalized vector of 

[n(s)J k {jf > 3 > is an eigenvector of the problem since it has been shown to equal ^cp( s)^> , 
provided that k is sufficiently large. By applying the orthogonality relationship and 
using equation (46a), the premultiplication by [n(s)J of the normalized eigenvector 
(cp(s)y obtained by the procedure just described yields 

[n(s)J {V( sj) = X(s) (cp(sj) (El 3) 

and thus, the corresponding eigenvalue is obtained. 

After each premultiplication of an arbitrary vector by ^N(s)jJ , the resulting vector 
may be normalized without affecting the convergence of the process. In the computational 
solution of the problem, the normalization of the vector after each premultiplication is an 
expedient method in that it permits a check on the convergence after each premultiplica- 
tion. From equation (El 3) it is known that a solution is obtained when the premultiplica- 
tion of the normalized vector by [n(s)|| yields the same normalized vector within cer- 
tain specified computational accuracy limits. 

The iteration procedure described for obtaining the sth mode may be used to deter- 
mine the natural vibration characteristics of each mode. Thus, the solutions of the 
eigenvalue problem are given by the iteration of equation (El 3) since it is proved that 
such iteration results in convergence to the sth mode. 
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VIBRATIONS OF A NONROTATING CANTILEVERED BEAM 


The differential equation of motion for a nonrotating beam vibrating freely in one 
plane is 


d 


2 


dx 2 





mw = 


0 


(FI) 


If the beam is assumed to be cantilevered, the boundary conditions are 


w 


|x=x f 


= 0 


dw 

dx 


= 0 


x=: 


x 0 


(F2a) 

(F2b) 



= 0 


x=x 


n 


(F2c) 


d_ 

dx 



dV\ 

dx 2 / 


x=x 


n 


= 0 


(F2d) 


The differential equation may be written for each of a chosen set of equally spaced 
stations defined by x^ = Xq + ih where h is the spacing interval and i = 0, 1, 2, . . .,n. 
Thus equation (FI) determines a set of (n + 1) equations which may be combined into a 
matrix equation as 

- w 2 [nf|{w} = {0} (F3a) 

where 






(w 


n 




(F3b) 
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(EI ) 0 






m 


Is 



(F3c) 


(F3d) 


Equation (F3a) may now be integrated four times by using the integrating matrix as 
an operator. However, after the second integration it is necessary to premultiply the 
equation by the inverse of the stiffness matrix. The result after each of the successive 
integrations is given in the following equations: 




{w} - o> 2 [l] 2 [m] {w} + kj[l]{l} + k 2 (l) = {0} 


<F5) 


~(w} - " 2 HW t'THw + k i[ i ][ s ] TOO) 

+ “MCI"'® + k 3<X> - 0} < P6 > 

(w> - u, 2 [,]ts] + kJY] [s] 

+ + k 3 H(l} + k 4 (l> = {0} (F7) 
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where are constants of integration. Also, {Y} and {o} are (n + 1) element column 
matrices with all elements 1 and 0, respectively. The inverse of the stiffness matrix is 
given by 



For the purpose of evaluating the constants of integration - kj, the following matrix 
operators which contain (n + 1) elements are defined. 


[_ B oJ = L 1 0 0 • • • °J 

[_ B nJ = |_0 ... 0 0 lj 


(F9a) 

(F9b) 


Multiplying equation (F4) by jjBjJ and applying the boundary condition given by 
equation (F2d) yields 


k i = 


(F10) 


Multiplying equation (F5) by [^B n J , applying the boundary condition of equation (F2c), and 
substituting equation (F10) gives k2 as 


k 2 = *\b„J[i] 2 [m] {w} - ^ 2 LBnJ[l]0> L^JCOD 5 ] 0) 


(Fll) 


By multiplying equations (F6) and (F7) by |^BqJ applying the remaining boundary 
conditions, it is found that 


k 3 = k 4 = 0 


(F12) 


The following eigenvalue problem results by substituting equations (F10), (Fll), 

o 

and (F12) into equation (F7), simplifying, and dividing by w : 


-L<w} =[i] 2 [s] WH® 

00 - [0> KJ - W]H 


where 


(FI 3a) 


(FI 3b) 
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TABLE 1.- PHYSICAL PROPERTIES OF PROPELLER BLADE 


(a) U.S. Customary Units 
Jxq = ^ in.; h = 2 inTj 


Station 

number, 

i 

m ’ / 

lb-sec2/in^ 

Ely'y'j 

lb -in 3 

EI Z’Z’’ 

lb-in 3 

13, 

deg 

0 

1.026 X 10-3 

0.200 x 106 

63 x 10 6 

30.5 

1 

.696 

.110 

49 

25.2 

2 

.660 

.083 

46 

20.0 

3 

.608 

.058 

44 

14.8 

4 

.564 

.042 

43 

9.6 

5 

.535 

.031 

43 

4.7 

6 

.520 

.027 

44 

0 

7 

.506 

.026 

47 

-4.2 

8 

.498 

.025 

51 

-7.6 

9 

_ _ . 

.498 

.024 , 

56 

-10.0 


(b) SI Units 

£xq = 0.1524 m; h = 0.0508 m] 


Station 

number, 

i 

m, 

N-sec 3 /m 3 

EI y'y’ ’ 
N-m^ 

EI z t 2 t j 

N-m^ 

deg 

0 

7.07 

0.574 x 10 3 

181 X 10 3 

30.5 

1 

4.80 

.316 

141 

25.2 

2 

4.55 

.238 

132 

20.0 

3 

4.19 

.166 

126 

14.8 

4 

3.89 

.121 

123 

9.6 

5 

3.69 

.089 

123 

4.7 

6 

3.59 

.077 

126 

0 

7 

3.49 

.075 

135 

-4.2 

8 

3.43 

.072 

146 

-7.6 

9 

3.43 

.069 

161 

-10.0 


59 



TABLE 2.- NATURAL VIBRATION FREQUENCIES OF PROPELLER BLADE 


/3 at 0.75R, 

Mode 

n, 

j CO, Hz 

Percent 

deg 

rpm 

Experimental 

Computed 

error 

0 

1 

1567 

40.08 

40.77 

1.7 



2609 

58.73 

59.85 

1.9 



3583 

76.52 

78.08 

2.0 



4486 

93.07 

95.03 

2.1 



5884 

117.50 

121.24 

3.2 


2 

1589 

107.53 

109.05 

1.4 



2614 

137.02 

139.03 

1.5 



3585 

170.60 

171.41 

.5 



4537 

202.53 

204.40 

.9 

20 

1 

1568 

37.66 

38.21 

1.5 



2585 

53.44 

54.32 

1.6 



4513 

83.74 

85.37 

1.9 



6016 

106.66 

109.46 

2.6 



6983 

122.08 

124.86 

2.3 


2 

1576 

107.42 

108.62 

1.1 



2561 

135.45 

137.26 

1.3 

40 

1 

1506 

33.10 

34.30 

3.6 



2530 

46.95 

47.65 

1.5 



4503 

72.28 

73.98 

2.4 



6007 

92.93 

94.01 

1.2 


2 

1484 

107.58 

106.05 

-1.4 



2532 

133.93 

135.71 

1.3 



3583 

167.69 

170.48 

1.7 


0 

0 


TABLE 3.- COMPARISON OF COMPUTED FIRST-MODE DISPLACEMENTS FOR 


n = 10, 25 AND r = 1, 2, . . .,7 WITH EXACT -SOLUTION VALUES 


x/Z 

Exact 

Number 

of 

intervals 



Degree of integrating matrix 



solution 

Seventh 

Sixth 

Fifth 

Fourth 

Third 

Second 

First 

0 

0 

25, 10 

0 

0 

0 

0 

0 

0 

0 

0.2 

0.12774 

25 

0.12774 

0.12774 

0.12774 

0.12774 

0.12774 

0.12774 

0.12754 



10 

.12774 

.12774 

.12774 

.12774 

.12775 

.12778 

.12649 

0.4 

0.45977 

25 

0.45977 

0.45977 

0.45977 

0.45977 

0.45977 

0.45978 

0.45946 



10 

.45977 

.45977 

.45977 

.45977 

.45977 

.45986 

.45782 

0.6 

0.92227 

25 

0.92227 

0.92227 

0.92227 

0.92227 ! 

0.92227 

0.92228 

0.92193 



10 

.92227 

.92227 

.92227 | 

.92227 

.92227 

.92237 

.92013 

0.8 

1.45096 

25 

1.45096 

1.45096 

1.45096 

1.45096 

1.45096 

1.45096 

1.45070 


L 

10 

1.45096 

1.45096 

1.45096 

1.45096 

1.45095 

1.45099 

1.44934 

1.0 

2.00000 

25, 10 

2.00000 

2.00000 

2.00000 

2.00000 

2.00000 

2.00000 

2.00000 


TABLE 4.- COMPARISON OF COMPUTED SECOND-MODE DISPLACEMENTS FOR 
n = 10, 25 AND r = 1, 2, . . ., 7 WITH EXACT-SOLUTION VALUES 


x/Z 

Exact 

Number 

of 

intervals 



Degree of integrating matrix 



solution 

Seventh 

Sixth 

Fifth 

Fourth 

Third 

Second 

First 

0 

0 

25, 10 

0 

0 

0 

0 

0 

0 

0 

0.2 

0.60211 

25 

0.60211 

0.60211 

0.60211 

0.60211 

0.60212 

0.60267 

0.60004 



10 

.60211 

.60210 

.60211 

.60255 

.60262 

.60926 

.58872 

0.4 

1.36694 

25 

1.36694 

1.36694 

1.36694 

1.36695 

1.36694 

1.36786 

1.36583 



10 

1.36693 

1.36694 

1.36699 

1.36757 

1.36667 

1.37777 

1.35990 

0.6 

1.17895 

25 

1.17895 

1.17895 

1.17895 

1.17896 

1.17895 

1.17944 

1.17933 



10 

1.17894 

1.17895 

1.17901 

1.17935 

1.17873 

1.18344 

1.18158 

0.8 

-0.14007 

25 

-0.14007 

-0.14007 

-0.14007 

-0.14007 

-0.14006 

-0.14010 

-0.13855 



10 

-.14007 

-.14008 

-.14004 

-.13990 

-.13925 

-.14047 

-.13003 

1.0 

2.00000 

25, 10 

2.00000 

2.00000 

2.00000 

2.00000 

2.00000 

2.00000 

2.00000 
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TABLE 5.- COMPARISON OF COMPUTED THIRD-MODE DISPLACEMENTS FOR 


n = 10, 25 AND r = 1, 2, . . .,7 WITH EXACT -SOLUTION VALUES 


x/l 

Exact 

solution 

Number 

of 

intervals 



Degree of integrating matrix 



Seventh 

Sixth 

Fifth 

Fourth 

Third 

Second 

First 

0 

0 

25, 10 

0 

0 

0 

0 

0 

0 

0 

0.2 

1.20901 

25 

1.20901 

1.20901 

1.20901 

1.20919 

1.20905 

1.21780 

1.20348 



10 

1.20819 

1.21004 

1.21164 

1.22304 

1.20864 

1.32008 

1.17216 

0.4 

1.05185 

25 

1.05185 

1.05185 

1.05185 

1.05196 

1.05180 

1.05735 

1.04898 



10 

1.05142 

1.05307 

1.05309 

1.05744 

1.04461 

1.11491 

1.03436 

0.6 

-0.94753 

25 

-0.94753 

-0.94753 

-0.94753 

-0.94760 

-0.94748 

-0.95061 

-0.95048 



10 

-.94717 

-.94842 

-.94905 

-.94985 

-.94014 

-.97675 

-.96540 

0.8 

-0.78975 

25 

-0.78975 

-0.78975 

-0.78975 

-0.78979 

-0.78981 

-0.79066 

-0.79543 



10 

-.78907 

-.79012 

-.79264 

-.79192 

-.79079 

-.78789 

-.82736 

1.0 

2.00000 

25, 10 

2.00000 

2.00000 

2.00000 

2.00000 

2.00000 

2.00000 

2.00000 


TABLE 6.- COMPARISON OF COMPUTED FOURTH-MODE DISPLACEMENTS FOR 
n = 10, 25 AND r = 1, 2, . . .,7 WITH EXACT -SOLUTION VALUES 


x/l 

Exact 

Number 

of 

intervals 



Degree of integrating matrix 



solution 

Seventh 

Sixth 

Fifth 

Fourth 

Third 

Second 

First 

0 

0 

25, 10 

0 

0 

0 

0 

0 

0 

0 

0.2 

1.50758 

25 

1.50757 

1.50762 

1.50761 

1.50914 

1.50737 

1.54968 

1.49965 



10 

1.51359 

1.54398 

1.52492 

1.58845 

1.46467 

2.04705 

1.45443 

0.4 

-0.63112 

25 

-0.63112 

-0.63113 

-0.63114 

-0.63160 

-0.63116 

-0.64391 

-0.63254 



10 

-.63215 

-.64456 

-.64410 

-.65788 

-.61936 

-.77565 

-.63891 

0.6 

-0.65299 

25 

-0.65299 

-0.65299 

-0.65301 

-0.65331 

-0.65303 

-0.66091 

-0.65156 



10 

-.65380 

-.66362 

-.66570 

-.66231 

-.64129 

-.72731 

-.64517 

0.8 

1.28608 

25 

1.28607 

1.28607 

1.28612 

1.28641 

1.28592 

1.29194 

1.29399 



10 

1.29344 

1.30838 

1.30332 

1.27247 

1.24597 

1.28186 

1.33918 

1.0 

2.00000 

25, 10 

2.00000 

2.00000 

2.00000 

2.00000 

2.00000 

2.00000 

2.00000 
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Natural vibration frequency, Ha 





Modal displacement 



Modal displacement, v Modal displacement, v 

(a) 0 rpm. (b) 6016 rpm. 

Figure 2.- Modal displacements at the equally spaced stations for the propeller blade having a pitch angle 0 

of 20° at x = 0.75R. 


Percent error 
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Figure 3.- Percent error of computed natural frequencies of uniform beam 
as a function of n for the first-degree integrating matrix. 
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Figure 4.- Percent error of computed natural frequencies of uniform beam 
as a function of n for the second-degree integrating matrix. 
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Figure 5.- Percent error of computed natural frequencies of uniform beam 
as a function of n for the third-degree integrating matrix. 
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Figure 8.- Percent error of computed natural frequencies of uniform beam 
as a function of n for the sixth-degree integrating matrix. 
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Figure 9.- Percent error of computed natural frequencies of uniform beam 
as a function of n for the seventh-degree integrating matrix. 
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